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ABSTRACT 

High speed photometric observations of the spectroscopically-discovered PG 1159 star 
SDSS J034917. 41-005917. 9 in 2007 and 2009 reveal a suite of pulsation frequencies 
in the range of 1038 - 3323 /iHz with amplitudes between 3.5 and 18.6 mmag. 
SDSS J034917. 41-005917. 9 is therefore a member of the GW Vir class of pulsating 
pre-white dwarfs. We have identified 10 independent pulsation frequencies that can 
be fitted by an asymptotic model with a contant period spacing of 23. 61 ± 0.21 s, 
presumably associated with a sequence of I = 1 modes. The highest amplitude peak 
in the suite of frequencies shows evidence for a triplet structure, with a frequency 
separation of 14.4 /iHz. Five of the identified frequencies do not fit the i — 1 sequence, 
but are, however, well-modeled by an independent asymptotic sequence with a con- 
stant period spacing of 11.66 ± 0.13 s. It is unclear to which i mode these frequencies 
belong. 
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1 INTRODUCTION 

The GW Vir stars are a sub-group of variables in the spec- 
troscopic PG 1159 class, which form a link between the 
(post-AGB) central stars of planetary nebulae and the H- 
L" ' deficient white dwarf cooling sequence. They pulsate non- 
. ,— | ■ radially and lie in an instability strip bounded by effec- 
^ ■ tive temperatures 200 000 > T cfI > 75 000 K, excited by 
S_h ' the kappa mechanism working through partial ionization of 
carbon and oxygen. Studying these stars with astroseismol- 
ogy has provided important knowledge on the interiors of 
the late stages of stellar evolution (Winget & Kepler 2008). 
There are 19 known GW Vir stars (Quirion, Fontaine & 
Brassard 2007; Quirion 2009), showing a wide variety of be- 
haviour. The possible addition of more examples is therefore 
of significance. Here we show that the known spectroscopic 
PG 1159 star SDSS J034917.41-005919.2 (hereafter SDSS 
J0349-0059) is a non-radial pulsator, putting it in the GW 
Vir subclass. In Sect. 2 we describe what is already known 
about SDSS J0349-0059 and list our high speed photometric 
observations. Sect. 3 analyses these and presents compar- 
isons with other GW Vir stars. 
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2 SDSS J0349-0059 AS A PG 1159 STAR 

2.1 The PG 1159 Instability Strip 

Hiigelmeyer et al. (2006) used Sloan Digital Sky Survey 
spectra of five spectroscopically discovered PG 1159 stars, 
including SDSS J0349-0059, to interpret their spectra us- 
ing non-LTE model atmospheres; these were added to the 
analyses of six previously analysed stars. Among these, 
SDSS J0349-0059 has an average position, with T cff = 90.0 
± 0.9 kK and logg = 7.50 ± 0.01 (cgs). 

Corsico et al. (2006) computed nonadiabatic pulsation 
models for stars in the GW Vir instability strip and found 
that their models agree with the observed strip. The param- 
eters listed above for SDSS J0349-0059 place it within the 
observed instability strip (cf. figures 6 and 7 of Corsico et 
al.). 



2.2 High speed photometry of SDSS J0349-0059 

SDSS J0349-0059 was observed in January 2007, March 2009 
and December 2009 with the University of Cape Town's 
CCD photometer (O'Donoghue 1995) attached to the 40- 
in and 74-in reflectors at the Sutherland site of the South 
African Astronomical Observatory. The motivation for ob- 
serving this star was inclusion as a UV-rich star included in 
a search for possible AM CVn stars. 

The observing log is given in Tab. 1 and the individual 
light curves are shown in Fig. 1. Note that the short ob- 
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Table 1. Observing log of photometric observations 



SDSS J0345917. 41-005919. 2 



20 



Run 

nun 


Date of obs. 


TT Tpit 

nj 1 ) 


Length 




v 

V 


CO 

e 




(start of night) 




W 


( s ) 


(mag) 


G 

c 


S7699 


20/01/2007 


4121.29287 


2.38 


20 


17.7 


-v 


S7701 


21/01/2007 


4122.27472 


3.39 


20 


17.7 


a 


S7703 


22/01/2007 


4123.27466 


3.02 


20 


17.7 


B 
< 


S7706 


23/01/2007 


4124.27198 


3.06 


20 


17.8 




S7846 


20/03/2009 


4911.23629 


1.34 


30 


17.9 




S7855 


23/03/2009 


4914.23355 


1.19 


30 


17.9 


CO 


S7859 


24/03/2009 


4915.23460 


1.12 


30 


17.9 


B 


S7899 


24/12/2009 


5190.35493 


1.27 


30 


17.8 




S7901 


25/12/2009 


5191.28264 


2.83 


30 


17.8 


XI 


S7904 


26/12/2009 


5192.28631 


2.63 


30 


17.8 


a 


S7906 


27/12/2009 


5193.28308 


2.74 


30 


17.9 


B 
< 



^HJD of first observation; tj n is the integration time. 



IB 



SDSS J034917. 41-005919. 2 

S7699 ^V^WjSJW^ 
S7701 ( + 0.3) tj*^^ 1 ^^*?*^*?^^ 
S7703 ( + 0.6) i^f^A^^f^^^^fi} 
S7706 ( + 0.9) ^A^^\Y^V^M^ 
\VV.- S7846 (+1.0) 



S7899 (+1.8) 



S7855 (+1.3) 
""^V S7859 (+1.6) 

S7901 ( + 2.1) »VV^^"'A i ''^H w V''WVv 
S7904 ( + 2.4) ^V^»VV , ^y - *N^V^^^v'<ys!«r/V 

S7906 ( + 2.7) *^l , Ja*^\\'-^A!'-r^*&^y?S 



0.25 



0.3 



0.35 



0.4 



Day 



Figure 1. Individual light curves of SDSS J0349-0059. The light 
curve of run S7699 is displayed at the correct brightness. Vertical 
offsets (indicated between brackets) have been applied to other 
light curves for display purposes only. 



serving runs in March 2009 (S7846, S7855 and S7859) were 
observed at high airmass (1.7 — 3); this is evident in the 
downwards trend seen in each of the light curves. The dif- 
ferential photometric correction fails at these large airmasses 
due to the different colours of the target and the reference 
star; SDSS J0349-0059 has (g-r) = -0.42, where the refer- 
ence star is redder at (g — r) — 0.26 (SDSS Data Release 8: 
Aihara et al. 2011). 

The Fourier transforms (FTs) of the three data sets 
(January 2007, March 2009 and December 2009) are shown 
in Fig. 2. In the construction of the FTs, individual runs 
have their mean and trend removed. 

Rapid oscillations with a range ~ 0.2 mag and time 
scale ~ 300 — 900 s are visible in the light curves and their 
FTs, which, with the spectra discussed by Hiigelmeyer et 
al. (2006), place SDSS J0349-0059 in the GW Vir subgroup 
of PG 1159 stars. 
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Figure 2. The Fourier transforms of SDSS J0349-0059 in January 
2007 (upper panel), March 2009 (middle panel) and December 
2009 (lower panel). 



3 ANALYSIS 

The Fourier transforms of the combined runs in January 
2007 (upper panel of Fig. 2) and December 2009 (lower panel 
of Fig. 2) show a substantial amount of variation in the 
distribution of power in the FTs. The March 2009 data set 
is shown for completeness, but given the relatively few short 
observing runs, the signal-to-noise in the FT of the combined 
runs in March 2009 (middle panel of Fig. 2) is too low to 
identify many individual pulsation frequencies. However, the 
main peak in March 2009 at ~ 2865 /iHz is also seen in 
January 2007 and December 2009. 

The distribution of pulsation frequencies in SDSS 
J0349-0059 resembles closely that of the GW Vir star 
PG 1707+427 (Kawaler et al. 2004), which has a main peak 
at 2236 /iHz and additional frequencies in the range of 1100 
- 2988 /kHz. In SDSS J0349-0059 the main peak is located at 
2386 /iHz with additional pulsation frequencies in the range 
of 1097 - 3323 /iHz. 

Fig. 3 shows an expanded view of the FT of the com- 
bined runs in January 2007 (upper two panels) and Decem- 
ber 2009 (lower two panels) in the 1000 - 3500 /iHz frequency 
range. The windows functions for the two data sets are given 
in the rightmost panels of Fig. 3, scaled to the amplitude of 
the highest peak in the respective FTs. In each FT, the hor- 
izontal dashed line corresponds to 4cr amplitude detection 
limit; this limit is 3.5 mmag for January 2007 and 3.1 mmag 
for December 2009. The standard deviation (a) has been 
determined over the 1000 - 4000 /iHz frequency range, after 
prewhitenening the FT by the frequencies listed in Tab. 2 
and randomizing the residual brightness variations for each 
observing run against the timing array. 

In the January 2007 FT, 13 frequencies have been iden- 
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Figure 3. The Fourier transforms of SDSS J0349-0059 in January 2007 (upper panels) and December 2009 (lower panels). For each data 
set three panels are displayed: the original Fourier transform of the combined observing run (top- left), the residual Fourier transform after 
prewhitening with the identified frequencies listed in Tab. 2 (lower-left) and the window function (top right). The 4a amplitude limit 
is indicated by the horizontal dashed line, the identified frequencies are marked by the vertical dashed lines and frequencies associated 
with the two asymptotic models with a constant period spacing are marked by short vertical markers, see Sect. 3.2 for details. 



tified with amplitudes above the 4er amplitude limit. They 
are marked by the vertical dash lines in Fig. 3. Successive 
prewhitening of the highest amplitude peak has resulted in 
the identification of the 13 pulsation frequencies listed in 
Tab. 2. Some of these frequencies could possibly be mis- 
identified as a one-day alias, although the coincidence of 
frequencies in the two independent data sets gives some re- 
assurance of correct identification; in Tab. 2 we merely list 
the frequencies associated with the highest amplitude sig- 
nals. The final FT prewhitened at these 13 frequencies is 
shown immediately below the original FT. 



Similarly, ten independent frequencies have been iden- 
tified in the December 2009 data. They are listed in Tab. 2 
and marked by the vertical dashed lines in the third vertical 
panel (left) of Fig. 3. Again, the final FT prewhitened at the 
identified frequencies is displayed below the original FT. 



The highest amplitude peak in both the 2007 and 2009 
observations is located at ~ 2386 /iHz (peak amplitude of 
18.6 and 15.2 mmag in 2007 and 2009, respectively). Other 
pulsation frequencies in common between the two data sets 
are around 2865 /iHz - peak amplitudes of 12.1 mmag (2007) 
and 9.6 mmag (2009) - and 3323 /iHz - peak amplitudes of 
5.5 mmag (2007) and 4.3 mmag (2009). 



3.1 Rotational splitting of the 2387-/iHz (419 s) 
oscillation 

The highest peak in the 2007 data shows clear evidence for 
three closely spaced frequencies. This is illustrated in Fig. 4 
where we show the January 2007 FT in the frequency range 
of 2000 - 3000 /iHz (original, top panel). After prewhitening 
by the main peak at 2387.2 /iHz, a clear window function re- 
mains (middle panel) at 2401.4 /iHz. Prewhitening at 2401.4 
/iHz leaves a distinct signal at 2372.6 /iHz (lower panel). All 
frequencies are marked by vertical dashed lines. This sug- 
gests that the strongest peak is in fact a triplet, with a 
frequency separation of 14.4 ± 0.2 /iHz. 

The 2009 data reveal two components of the triplet; a 
strong peak at 2385.6 /iHz and a low amplitude peak at 
2383.5 /iHz which is probably the 1-day alias of the 2372.6 
/iHz peak seen in 2007. 

We draw again on the similarity with PG 1707+427, 
where evidence was found for a 9 /iHz split of the main 
mode, interpreted as a rotational frequency splitting, cor- 
responding to a rotation period of 0.65 days (Kawaler et 
al. 2004). In SDSS J0349-0059 this would imply a rotation 
period of 0.40 ± 0.01 days if the pulsation mode is £ = 1. 

3.2 Period spacing 

Given the strong similarity with PG 1707+427, we suspected 
that the frequencies identified in SDSS J0349-0059 and listed 
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Table 2. Frequency identifications 
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Amplitude Rank 



Table 3. Periods compared to an asymptotic model with a period 
spacing of 23.61 s and a constant period of 300.54 s. 
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Figure 4. The Fourier transform of SDSS J0349-0059 in Jan- 
uary 2007 (upper panel). The middle panel displays the Fourier 
transform after prewhitening of the strongest signal in the orig- 
inal data at 2387.2 /iHz. The lower panel displays the residual 
Fourier transform after subsequently prewhitening the data at a 
frequency of 2401.4 /tHz. Peak frequencies of the resolved triplet 
at 419 s are marked by the dashed vertical lines. 



in Tab. 2 are predominantly part of a sequence of I = 1 g- 
modes with constant period spacing; this is an asymptotic 
model of the form P mo d e i = nAP + C, where P mo dei is the 
model period, n is an integer, AP is a constant period spac- 
ing and C is a constant period. 

Guided by the analysis of Kawaler et al. (2004), we first 
looked at the pulsation frequencies in the December 2009 
data. In the December 2009 observations, a large-amplitude 
low frequency pulsation was detected at 1097 /iHz, assumed 



3 The year in which the pulsation is detected; the amplitude of 
the pulsation in listed in square brackets (in units of mmag). 

* The m = ±1 modes of this triplet are not listed here. 
*The one-day alias of the An = 8 mode is at 489.2 s, which 

gives 8P = -0.2 s and 5P/23.61 = -0.01. 
*The one-day alias of the An = 26 mode (2007) is at 916.2 s, 
which gives SP = 1.8 s and 5P/23.61 = 0.07. 



to be a high order (n) mode similar to the An = 25 mode 
identified by Kawaler et al. in PG 1707+427 around 1100 
/iHz. 

To search the suite of periods for a constant period 
spacing, we used the inverse variance method as outlined 
by O'Donoghue (1994). By letting AP vary between 10 and 
40 s (in steps of 0.005 s), we compared the measured periods 
with the model periods and plotted the inverse variance of 
the timing residual (divided by AP) as a function of AP. 
The results are shown in the upper panel of Fig. 5 for all 
the 10 periods identified in December 2009. A low ampli- 
tude peak around AP ~ 23.5 is evident. The significance 
of this peak improves substantially when two periods are 
removed from the sequence (412.3 and 516.7 s), where the 
remaining 7 independent periods (we have also removed the 
lower sideband component of the 419-s triplet) form a well- 
defined sequence with AP = 23.52 ± 0.19 s and C = 301.70 
s. This is shown in the bottom panel of Fig. 5. In compari- 
son, Kawaler et al. (2004) find AP = 23.0 s and C = 332.9 
s for PG 1707+427. 

All but five periods identified in 2007 and 2009 fit the 
asymptotic model with a constant period spacing around 
~23.5 s. They are listed in Tab. 3 against the identified mode 
difference (An), the period difference (SP = P b s — Pmodei) 
and presumably belong to £ = 1 modes. The best fit model 
based on the periods listed in Tab. 3 gives AP = 23.61 + 0.21 
s and C = 300.54 s. The upper panel of Fig. 6 shows the 
inverse variance as a function of AP for this set of periods. 
Had we only selected frequencies with amplitudes above a 
5(J limit (as opposed to our current 4a selection), the inverse 
variance analysis would have given a similar result (AP = 
23.47+0.13 s and C = 301.67 s). The sequence of frequencies 
associated with the best fit model (AP = 23.61 s) is shown 
in Fig. 3 by the lower sequence of short vertical markers 
(starting at 3327 /iHz). 

The remaining five periods not matched by this model 
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Figure 5. The inverse variance versus AP for all periods identi- 
fied in December 2009 (upper panel), and for seven independent 
periods (all except 412.3, 419.6 and 516.7 s; lower panel). 



Table 4. Periods compared to an asymptotic model with a period 
spacing of 11.66 s and a constant period of 353.91 s. 
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§ The year in which the pulsation is detected; the amplitude of 
the pulsation in listed in square brackets (in units of mmag). 

are part of a sequence themselves with a fixed period spac- 
ing of AP = 11.66 ± 0.13 s and C = 353.91 s. They arc 
listed in Tab. 4 and the result of the inverse variance test 
is shown in the lower panel of Fig. 6. Within the formal 
error, this period spacing is consistent with being equal to 
half the period spacing of the 1 = 1 model. The sequence 
of frequencies associated with this model (AP = 11.66 s) is 
also shown in Fig. 3 by the upper sequence of short vertical 
markers (starting at 2826 /kHz) and can be compared with 
the asymptotic model for the £ = 1 modes. 

The 8P/AP residuals we obtain for the two models are 
similar to those found in studies of other GW Vir stars (e.g., 
Kawaler et al. 2004). 



4 DISCUSSION 

Our photometric observations show that SDSS J0349-0059 is 
a non-radial pulsator in the GW Vir class of variable stars. 
Frequency splitting of the principal oscillation mode at 419 
s reveals a rotation period of 0.40 d, if the oscillation mode 
is I = 1. As with many other GW Vir stars, it is possible 
to represent most of the modes with a linear relationship in 
period, similar to what has been seen in PG 1707+427, but 
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Figure 6. The inverse variance versus AP for all the periods 
identified as £ = 1 modes (upper panel) and for the remaining set 
of 5 periods that cannot be fitted by the £ = 1 constant spacing 
model (lower panel). 



with parameters that put SDSS J0349 slightly redward of 
the latter star, which now defines the red edge of the GW 
Vir instability strip. 

The five oscillation modes not included in the above 
group can be fitted by another linear relationship in which 
the spacing is a harmonic of the first sequence. However, 
there is no currently known physical model which explains 
this behaviour. This is unlike PG 1707+427 where the dis- 
crepant oscillation modes appear to have an £ — 2 origin. 
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